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We study the g-state Potts model on the simple cubic lattice with ferromagnetic interactions in one 
lattice direction, and antiferromagnetic interactions in the two other directions. As the temperature 
T decreases, the system undergoes a second-order phase transition that fits in the universality class 
of the 3D 0(n) model with n — q — 1. This conclusion is based on the estimated critical exponents, 
and histograms of the order parameter. At even smaller T we find, for g = 4 and 5, a first-order 
transition to a phase with a different type of long-range order. This long-range order dissolves at 
T = 0, and the system effectively reduces to a disordered two-dimensional Potts antiferromagnet. 
These results are obtained by means of Monte Carlo simulations and finite-size scaling. 
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I. INTRODUCTION 

According to the hypothesis of universality, critical 
phase transitions fall in classes determined by spatial di¬ 
mensionality and symmetry of the order parameter. The 
latter is usually reflected by the degeneracy of the ground 
state of the Hamiltonian. However, for certain systems 
at criticality, a higher symmetry may emerge in the or¬ 
der parameter, and the associated critical behavior may 
become very rich. 

Examples of this phenomenon are known, in particu¬ 
lar in two dimensions (2D). In the g-state clock modeli, 
the spins are confined to a plane and take q discrete di¬ 
rections (j) = 2Trnlq. The ferromagnetic ground state 
is g-fold degenerate, reflecting the Zq symmetry. For 
g > 5, however, the 2D clock model exhibits a Berezinskii- 
Kosterlitz-Thouless (BKT) transition as the temperature 
T is lowered, and quasi-long-range order with continu¬ 
ous 0(2) symmetry emerges at T < Tc, just as in the 
rotationally invariant XY model. Emergent symmetries 
are also found in many other physical systems, including 
a spin ice system^, deconfined quantum critical points^, 
high-Tc superconductors^ and so forth^, and are often ac¬ 
companied by very interesting phenomena. For example, 
in the spin ice^, emergent SU(2) symmetry leads to an 
unusual phase transition with a jump in the order pa¬ 
rameter, which is a feature of discontinuous transitions, 
whereas the the domain wall tension vanishes, which is 
a feature of continuous transitions. In a class of models 
with Z 2 and U(l) symmetry^, emergent supersymmetry 
at the Ising-BKT multicritical point leads to new critical 
behavior, with unusual scaling of the correlation length. 

The Potts model^ has spins with q values a = 
1,2,-- - ,q that interact as K6a-i,a-j, reflecting permuta¬ 
tion symmetry. The antiferromagnetic g = 3 model on 
the simple cubic lattice breaks an effective Zq symmetry 
at low temperatures, and 0(2) symmetry emerges at the 
critical point^i^. In line with these findings, a BKT tran¬ 
sitions with emergent 0(2) symmetry can arise on simple 
cubic lattices with a finite thickness^. For similar models 
with g > 3 one might thus expect emergent 0(n) symme¬ 
try, i.e., isotropy in n dimensions. While this scenario is 
consistent with numerical resultsiS, the situation is not 


entirely clea r^°d^ . 

It has been hypothesized^ that in high-Tc superconduc¬ 
tors, the magnetic and superconducting degrees of free¬ 
dom can merge into a critical state with effective SO(5) 
symmetry. However, as argued by Fradkin et al^, the 
symmetry of the corresponding 0(5) fixed point is eas¬ 
ily broken since the components of the order parameter 
are inequivalent on the microscale. This applies as well 
to other systemsi^ for which higher emergent symmetries 
had been proposed. 

A different critical behavior occurs in two-dimensional 
systems with mixed interactions—i.e., ferromagnetic 
(FM) in one direction and antiferromagnetic (AF) in the 
other. The g = 3 mixed Potts model on the square lat¬ 
tice undergoes a BKT-like transition, and 0(2) symmetry 
emerges in the low-temperature rangei^. Rich phenom¬ 
ena also occur in the the mixed Ising (g = 2) model on 
the multi-layered triangular latticei^. 

In this paper, we study g-state mixed Potts models on 
the simple cubic lattice, with FM couplings in the 2 : direc¬ 
tion and AF couplings in the x — y plane. Using cluster- 
type Monte Carlo algorithms, we find continuous phase 
transitions for 2 < g < 6 as the temperature T is lowered. 
The critical behavior of these systems is consistent with 
0(n) universality in 3D, with n = g — 1. This result may 
hold more generally, i.e., for at least some g > 6. To our 
knowledge, such an emergent 0(n) symmetry in the 3D 
Potts model has, apart from the g = 3 antiferromagnet, 
not been reported in literature. 

II. MODEL, ALGORITHM, CRITICAL POINTS, 
AND CRITICAL EXPONENTS 

The reduced Hamiltonian of the mixed Potts model is 

^ ^E E ^E E> (1) 

z <i,i> 2 : I 

where the spins take g values cr = 1,2, ••• ,g. With 
K > 0, the sum in the first term, taken over all nearest- 
neighbor sites {i,j) in layer z, defines AF couplings. The 
second term defines FM couplings in the z direction. We 
refer to the temperature as T = 1/K. 
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Cluster Monte Carlo methods are very effective for sim¬ 
ulation of FM lattice modelsi^, while the efficiency of 
the Wang-Swendsen-Kotecky (WSK) algorithm^ for AF 
Potts models depends on the lattice type and tempera¬ 
ture. For mixed interactions, we apply a single-cluster 
algorithm merging elements of the Wolff method for FM 
modelsi^) and the WSK algorithm. A combination with 
the geometric cluster algorithntii^— , which employs lat¬ 
tice symmetries, is still needed for effective simulations 
of systems up to L = 128 at sufficiently low tempera¬ 
tures. In addition we applied Metropolis sweeps. 

The sampled observables include the staggered suscep¬ 


tibility Xs, the uniform susceptibility Xu, their Binder 
ratios Qu and Qs, and the specific heat C^'- 

Xs = V{M!), Q. = (2) 

Xu = V{Ml), Qu = (3) 

Cv = V{{£^) - {£)^)/T\ (4) 

where V = is the system volume, and £ the energy 
density. and are defined as 

- Pb,pf (5) 

^ p=i 

■^1 =-^J2^Pa,p +Pb,pf -^ ( 6 ) 

^ p=i ^ 


with (with k = a, b) the density of state-p spins on 
sublattice k, namely 


1 

Pk,p — y 


rGk 


(7) 


with spin coordinates r = {x,y,z). The sublattice k is 
defined by the parity oi x + y. The Ms and Mu are 
the order parameters of the model, exposing a possible 
symmetry breaking of the model. It should be noted 
that, in the AF Potts model, the type of order in the 
low-temperature phases depends essentially on entropy 
effects, apart from the energy effect. For example, the 
low temperature phase of the 3-state AF Potts model on 
the simple cubic lattice displays long-range order with 
one sublattice frozen in one of the Potts states, while the 
spins on the other sublattices are free to randomly take 
one of the other Potts states^. The maximal entropy of 
the latter sublattice explains the existence of this type of 
state. The phase transition to this state is thus, at least in 
part, entropy-driven, similar to behavior found for certain 
two-dimensional q = 4 Potts antiferromagnets^. 

We have investigated the model ([T]) for g = 3,4,5 
and 6, with periodic boundary conditions. The proce¬ 
dure involved three steps, specified here for q = A. 
First, we simulate for several L at a number of tem¬ 
peratures T taken in a wide range. Each data point is 
based on 2 x 10® Monte Carlo steps (MCSs). Each MCS 
consists of 5 Wolff-cluster updates of the WSK type, 5 
geometric-cluster updates, 5 Metropolis sweeps, and data 
sampling. Each different simulation uses a different ran¬ 
dom seed, and starts from a random initial configuration, 
after which about 5x 10® data samples are discarded to al¬ 
low for equilibration of the system. These simulations are 


distributed over different CPU cores. After their comple¬ 
tion, the resulting data are collected, and the averages of 
the physical variables and their error bars are calculated. 
Plots of {Mg) and Qs in Fig.[T]yield an approximate crit¬ 
ical point Tc ~ 0.91. This is also seen in the scaling of 
the specific heat, in the left panel of Fig. [H 



T T 


FIG. 1. (Color online) Squared staggered magnetization den¬ 
sity {Ms) (left) and dimensionless ratio Qs (right) versus T 
for the 4-state mixed Potts model. The staggered magnetiza¬ 
tion becomes non-zero below Tc « 0.91. Near T « 0.55, {Ma) 
displays a jump, signaling a first-order phase transition. For 
T < 0.2, the correlation length in the z direction exceeds the 
system size, and crossover occurs to the disordered T = 0 
ground states. 



FIG. 2. (Color online) Specific heat CV (left) and squared 
magnetization density (Mu) (right) versus T for the 4-state 
mixed Potts model. The specific heat Cv displays a cusp at 
Tc « 0.91, while {M^) develops a discontinuity at the first- 
order transition ~ 0.55. 

Next, in order to determine the critical exponents of 
this phase transition we simulate near the estimated Tc, 
with 1.6 X 10^ MCS taken at each data point. The Qs 
data scale as 

where yt > 0 is the thermal exponent, ?/i < 0 is the 
correction-to-scaling exponent, and Qso, Ufe and b are un¬ 
knowns. A least-squares fit of this formula to the data 
yields Tc = 0.91381(1), and yt = 1.402(6). 

In the last step, we simulate at Tc and fit the data of 
Xs by 

Xs = L'^y''-^{a + bLy^), (9) 

with spatial dimensionality d = 3. The magnetic expo¬ 
nent follows as yh = 2.485(3). The uniform susceptibility 
Xu was also fitted by Eq. (0), with yh replaced by another 
magnetic exponent y^- This fit gives y^ = 1.787(5). 

The results for <7 = 3,5 and 6 used the same procedure. 
The critical points and exponents are listed in Table HI A 
comparison with the exponents yt and yu of the 0(n) 
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TABLE I. Critical points and critical exponents of the g-state mixed Potts model on the simple cubic lattice; a reflects a 
conservative average of some recent results given in Refs.[^ andl^. The critical exponents yt and yu of the recent results for 
the 3D 0(n) model, which are taken from Ref. [^, are also listed here for comparison. 


<? 

Tc 

g-state mixed Potts model 
yt Vh 

Vu 

n 

0(n) model 

yt 

yh 

2 

2.2557616(8)“ 

1.5873(1)“ 

2.481846(15)“ 

= yt 

1 

1.588(2) 

2.483(3) 

3 

1.36086(1) 

1.488(4) 

2.483(2) 

1.754(4) 

2 

1.488(3) 

2.483(2) 

4 

0.91381(1) 

1.402(6) 

2.485(3) 

1.787(5) 

3 

1.398(2) 

2.482(2) 

5 

0.64116(2) 

1.33(2) 

2.482(3) 

1.806(5) 

4 

1.332(7) 

2.483(2) 

6 

0.44545(2) 

1.30(3) 

2.485(3) 

1.837(4) 

5 

1.275(12) 

2.483(3) 


model2^ shows that the phase transition of the g-state 
mixed Potts model fits the universality class of the 0(n) 
model with n = q — 1. 


III. HISTOGRAM OF THE 0(N) SYMMETRY 

The remarkable emergence of 3D 0(n = g — 1) univer¬ 
sality in these models invites the construction of order 
parameter histograms, by representing the spins as vec¬ 
tors. These q vectors are symmetrically distributed in 
(g — l)-dimensional space, such that their scalar prod¬ 
uct matches the pair potential of Eq. which can be 
written as 


ecr • Ccr = 1, (no summationon a) (10) 

ea ■ Scr' =-(cr 7^ cr') (H) 

g- 1 

For example, in the case of g = 3 these vectors span a 
regular triangle, z.e., Co- = (cos0,sin0) with 9 = 2cr7r/3. 
For g = 4, the vectors are three-dimensional ones, and 
span a regular tetrahedron: 

Co-= (-l-l,-l-l,-l-l)/v^, for CT = 1; 

= (+1,-1,-1)/V3, fora = 2; 

= (-1,-1,-I-1)/a/ 3, for cr = 3; 

= (-1,-1-1,-1)/V^, for CT = 4.. 

Based on this vector representation, the magnetization is 
sampled separately for sublattices a and b. This yields 
the components of Ms and AIu in Eqs. m and (IH) 

Ms = Ma-Mb, Mn=Ma+Mb- ( 12 ) 

Figures I3i:a)-(c) display the histograms of the stag¬ 
gered magnetization for L = 32 systems, projected on 
two Cartesian axes, for g = 3, 4, and 5 respectively. The 
histograms are the same for any choice of the axes. The 
apparent isotropy shows the emergent 0(g — 1) symmetry. 
For g = 3, the symmetry persists in a finite system for 
a range below Tc, as shown in Fig. [3Kd). Figure |4] shows 
that Potts ferromagnets behave differently. It compares 
the histogram of the orientation (j) of the staggered mag¬ 
netization of the g = 5 mixed Potts model, projected on 
the a:, y plane, to similar plots for the magnetization of 
g = 3 ferromagnets with d = 2 and 3. 



FIG. 3. (Color online) Histograms of Ms for (a) g = 3 at 
criticality; (b) g = 4 at criticality; (c) g = 5 at criticality; and 
(d) g = 3 at T = 1.30 < Tc = 1.36. 



FIG. 4. (Color online) Histogram of the orientation cj) of Ms 
of a critical L = 32, g = 5 mixed Potts model (blue line). It 
is, modulo tt, given by </> = arctan(A4sy/A4sx). It is compared 
with analogous data for A4u of g = 3 Potts ferromagnets at 
Tc with d = 2 (T = 240, red peaks) and with d = 3 (L = 32, 
green peaks). The d = 3 transition is weakly discontinuous. 
Its data, shifted by 7r/3 for clarity, display a background from 
the coexisting disordered phase. The histograms are rescaled 
to match each other. 
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IV. FIRST ORDER TRANSITION OF THE 
MODEL 

In the range below Tc, the 4-state mixed Potts model 
displays jumps in Ms and A4u, near « 0.55 in Figs. [T] 
and[21 In the middle range <T < Tc, {M^ vanishes 
for L — 7 ^ oo, while {Ml) converges to a nonzero value. 
In contrast, both {Ml) and {M^ converge to nonzero 
values in the low T range 0 < T < T^. Thus differ¬ 
ent symmetries are broken on the two sides of T^,. The 
histograms of AIs and A4u show a broken Z 2 x 5'4 sym¬ 
metry for 0 < T < Tc, with Z 2 for the permutation of 
the two sublattices, and S 4 for the symmetric group for 
the four-state Potts model. In a typical configuration at 
0 < T < Tc, one Potts state dominates one sublattice, 
and the remaining states randomly occur on the other 
sublattice. This implies the breaking of the corner-cubic 
symmetry described by the four vectors for g = 4, pre¬ 
ceded by a sublattice sign ±. In a typical configuration 
at Tc < T < Tc, one sublattice is dominated by two ran¬ 
dom spin states, and the other sublattice by the other 
states. The staggered magnetization vector then points 
at one face of a cube, signaling a broken face-cubic sym¬ 
metry. The histograms of Ms are shown in Fig. [5ja) for 
0 < T < Tc and Fig. (S^b) for T^ < T < Tc, clearly dis¬ 
playing the broken corner-cubic and face-cubic symme¬ 
tries. Similarly, Figs. [S^d) and [5je) show the histogram 
of AIu for 0 < T < Tc, and T^ < T < Tc, respectively. 
The transition at T^ is not well visible in the ordinary en¬ 
ergy density E = {£), but it is exposed by the energy-like 
quantity Ei = {£{) based on the next-nearest neighbor 
correlations in the x-y planes, expressed as 

<<i,3» 

where << i,j » denote the next-nearest neighboring 
sites in the X — Y planes. Figure [S] shows the curves of 
E versus T and the curves of Ei versus T, the curves 
of E are featureless for this transition, but the curves of 
El at T'c for large system sizes show an obvious energy 
gap. This result reflects the stronger next-nearest neigh¬ 
bor correlations in the corner-cubic phase. 



q=4. L=4S, T=0.4 

(a) 


q=4. L=48, T=0.7S 

(b) 


q=3, L=J2, T=OJ 
(C) 




q=4, L=4S, T=0.4 

(d) 


q=4. L=4S, T=0.7S 

(e) 


q=3, L=32. T=OjS 

(f) 


FIG. 5. (Color online) Histogram of Ms and Mu for the 4- 
state and 3-state mixed Potts models. Lines are added to 
guide the eyes, especially to a 3D impression. The arrows 
describe the symmetry of Ads and of Mu- 



T 
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FIG. 6. (Color online) Left: energy density E of the 4-state 
mixed Potts model. Right: energy-like quantity density Ei, 
defined in the text, of the 4-state Potts model. The behav¬ 
ior of the curves of E\ with increasing system sizes clearly 
indicates the existence of a discontinuity at ~ 0.55 in the 
thermodynamic limit. 


The 3-state mixed Potts model breaks the Z 2 x S '3 
symmetry in the whole range 0 < T < Tc, as shown in 
Figs.[5Kc) and (f). But for the g = 5, another similar dis¬ 
continuous transition appears at T^ ~ 0.4, which breaks 
a Z 2 X S 5 symmetry in the low-T range 0 < T < T^- 
The degeneracy in the intermediate range T^ <T < Tc is 
2 X Cf = 20, where Cf denotes the binomial coefficient. 
In a typical configuration at T^ < T < Tc, the spins on 
one sublattice randomly take two states, and the other 
spins randomly take the remaining three states. 

A similar discontinuous transition may occur in the 6 - 
state mixed Potts model. We observed that the Z 2 x Sq 
symmetry is broken at low T, and that another ordered 
phase exists at intermediate T. But we did not find a 
jump in Ms or Mu for systems up to T = 96. This may 
still be due to a strong finite-size effect. 

At zero temperature the model reduces to a T = 0 
square-lattice AF Potts model^i, which is Neel ordered 
for g = 2, critical for g = 3, and disordered for g > 3. 
Figure [7] summarizes the phase behavior of the g-state 
mixed Potts models. 
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FIG. 7. (Golor online) Phase diagram of the g-state mixed 
Potts model. The multiplicities of the intermediate phases 
are expressed with help of the binomial coefficients Cl and 
Cf. The question mark means that the phase has not been 
unambiguously identihed. The models are disordered at T = 0 
for g > 3. 


V. CONCLUSION AND DISCUSSION 


In summary, our results indicate that g-state Potts 
models on the simple cubic lattice with mixed FM and AF 
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interactions display continuous phase transitions, with 
critical exponents in the 0(n) universality class with 
n = q — 1. The order parameter displays this emergent 
symmetry at criticality. In the low temperature ranges 
of the q = 4 and 5 models, perhaps also for g > 5, a dis¬ 
continuous transition occurs between two ordered phases. 
For T —>■ 0, the model crosses over to the d = 2 square- 
lattice AF Potts model, which is disordered for q > 3. 

Although the temperature is the only variable, the q- 
state mixed Potts model displays diverse and enigmatic 
phenomena. The 0(n) symmetry is not at all obvious in 
the Hamiltonian, but it nevertheless emerges, and con¬ 
trols the critical properties of the continuous phase tran¬ 
sition at Tc- In the sense of universality, the 3-state 
mixed Potts model is similar to the 0(2) model with a 
Zq perturbation^^, which also displays an emergent 0(2) 
symmetry at criticality. For q > 3, the mixed Potts model 
has a low-temperature ordered phase that spontaneously 
breaks the {q— I)-dimensional face-cubic symmetry. The 
histogram of the order parameter at criticality shows an 
emergent 0 (n = q—1) symmetry, and the estimated ther¬ 
mal exponent yt is a decreasing function of q, consistent 
with the 0(g—1) universality. In the analogous case of 
the pure antiferromagnetic g = 4 Potts model, the low- 
temperature ordered phase also breaks the 3D face-cubic 
symmetry and that the Monte Carlo simulation up to 
L = 96 also yields critical exponents consistent with the 
0(3) universality classic. Since the cubic perturbation is 


expected to be relevant for the 0(3) model^S, one cannot 
fully exclude that for g = 4, these phase transitions are of 
weak first order or belong to another universality class. 
However, in either case the question still remains why the 
effects of this perturbation are invisible in our analysis of 
finite systems. On the basis of our systematic study of 
the g-state mixed Potts models, we conclude that any 
symmetry-lowering perturbations of the emergent sym¬ 
metries are strongly suppressed, allowing the possibility 
that the ordering transitions fit exactly in the 0(g — I) 
universality classes. 

Finally, we mention that the mixed Potts model resem¬ 
bles a square-lattice quantum Potts antiferromagnet in a 
transverse fieldSIr— . The z dimension in the classical 
model corresponds with imaginary time in the Suzuki- 
Trotter formulation of the quantum model. The present 
series of mixed Potts models may provide a simple exam¬ 
ple where quantum fluctuations give rise to rich behavior. 
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